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In the framework of the isospin-dependent Boltzmann-Uehling-Uhlenbeck transport model, effect 
of the momentum dependence of nuclear symmetry potential on nuclear transverse and elliptic flows 
in the neutron-rich reaction 132 Sn+ 124 Sn at a beam energy of 400 MeV/nucleon is studied. We find 
that the momentum dependence of nuclear symmetry potential affects the rapidity distribution of 
the free neutron to proton ratio, the neutron and the proton transverse flows as a function of rapidity. 
The momentum dependence of nuclear symmetry potential affects the neutron-proton differential 
transverse flow more evidently than the difference of neutron and proton transverse flows as well as 
the difference of proton and neutron elliptic flows. It is thus better to probe the symmetry energy 
by using the difference of neutron and proton flows since the momentum dependence of nuclear 
symmetry potential is still an open question. And it is better to probe the momentum dependence 
of nuclear symmetry potential by using the neutron-proton differential transverse flow the rapidity 
distribution of the free neutron to proton ratio. 
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I. INTRODUCTION 

Nowadays the equation of state (EOS) of isospin sym- 
metric nuclear matter is now relatively well determined 
mainly by studying collective flows in heavy-ion collisions 
and nuclear giant monopole resonances [l|, [2| . The ma- 
jor remaining uncertainty about the EOS of symmetric 
nuclear matter is due to our poor knowledge about the 
density dependence of the nuclear symmetry energy [l|,|3|- 
|fj], which is crucial for understanding many interesting 
issues in both nuclear physics and astrophysics |7HlO|. 
And it is also crucial in connection with the structure 
of neutron stars and the dynamical evolution of proto- 
neutron stars [ll(. Nowadays considerable progress has 
been made recently in determining the density depen- 
dence of the nuclear symmetry energy around the nor- 
mal nuclear matter density. However, much more work is 
still needed to probe the high-density behavior of the nu- 
clear symmetry energy. Currently, to pin down the sym- 
metry energy, the National Superconducting Cyclotron 
Laboratory (NSCL) at Michigan State University, the 
Gesellschaft fuer Schwerioncnforschung (GSI) at Darm- 
stadt, the Rikagaku Kenkyusho (RIKEN, the Institute 
of Physical and Chemical Research) of Japan, and the 
Cooler Storage Ring (CSR) in Lanzhou are planning to 
do related experiments to probe the symmetry energy. 

The neutron-proton differential transverse flow and the 
difference of neutron and proton flows are both sensitive 
to the symmetry energy 0, [j], but the used transport 
models always adopt different momentum dependent in- 
teractions among nucleons. The importance of the mo- 
mentum dependence of nuclear symmetry potential on 
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the two kinds of nuclear flows was seldom mentioned. 
Considering the momentum dependence of nuclear sym- 
metry potential is quite controversial (12J . in the frame- 
work of the isospin-dependent Boltzmann-Uchling- Uh- 
lcnbcck trans por t model, we find that, besides the ra- 
tio of 7r~ /tt + [13[ , the momentum dependence of nuclear 
symmetry potential affects the neutron-proton differen- 
tial transverse flow more evidently than the difference of 
neutron and proton transverse flows as well as the differ- 
ence of proton and neutron elliptic flows. It is thus better 
to probe the symmetry energy by using the difference of 
nucleonic transverse (or elliptic) flows. 



II. THE IBUU04 TRANSPORT MODEL 

The present study is based on the IBUU04 transport 
model Q . The initial neutron and proton density distri- 
butions of the projectile and target are obtained by using 
the relativistic mean field theory. Although nuclear flow 
be affected by the initializations of colliding nuclei 
our main results in the present studies are not sen- 
sitive to the initializations of target and projectile nuclei. 
The experimental free-space nucleon-nucleon (NN) scat- 
tering cross sections and the in-medium NN cross sec- 
tions can be used optionally. In the present work, we 
did not use the isospin-dependent in-medium NN elas- 
tic cross sections from the scaling model according to 
nucleon effective masses [H, [lj| , although it is similar 
to the Brueckncr approach calculations p^7H19J . This 
is because the momentum-independent symmetry poten- 
tial (MID) can not use the isospin-dependent in-medium 
NN elastic cross sections from the scaling model accord- 
ing to nucleon effective masses. Thus all the studies 
here optionally use the experimental free-space nucleon- 
nucleon (NN) scattering cross sections for consistence. 



For the inelastic cross sections people use the experi- 
mental data from free space NN collisions since the in- 
medium inelastic NN cross sections are still very much 
controversial. The total and differential cross sections 
for all other particles are taken either from experimental 
data or obtained by using the detailed balance formula. 
In the model, besides nucleons, A and N* resonances 
as well as pions and their isospin-dependent dynamics 
arc included. The isospin dependent phase-space distri- 
bution functions of the particles involved are solved by 
using the test-particle method numerically. The isospin- 
dependence of Pauli blockings for fermions is also consid- 
ered. The momentum-dependent single nucleon potential 
(MDI) adopted here is [20J 
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In the above equation, 5 = (p„ — p p )/{p n + Pp) is the 
isospin asymmetry parameter, p = p n + p p is the baryon 
density and p n ,P P are the neutron and proton densi- 
ties, respectively, r = l/2(— 1/2) for neutron (proton) 
and r 7^ r' , a — 4/3, / r (r,p) is the phase-space dis- 
tribution function at coordinate r and momentum p. 
The parameters A u (x), Ai(x),B, C T)T , C TjT i and A were 
set by reproducing the momentum-dependent potential 
U(p, S, p, t) predicted by the Gogny Hartree-Fock and/or 
the Brueckner-Hartree-Fock calculations, the saturation 
properties of symmetric nuclear matter and the symme- 
try energy of about 32 MeV at normal nuclear matter 
density pq = 0.16 fm -3 . The propagations of nucleon are 
according to Hamilton's equations 



of momentum dependence of nuclear symmetry potential 
on the transverse and elliptic flows, we let the variable 
x be 1, since the IBUU04 model give s a super-soft sym- 
metry energy at higher densities [21| . In fact, behavior 
of nuclear symmetry energy at supra-densities is still in 
controversy. The main characteristic of the present single 
particle is the momentum dependence of nuclear symme- 
try potential, which has evident effect on energetic free 
n/p ratio in heavy-ion collisions [22| ■ In the present stud- 
ies, to show the effects of the momentum dependence of 
symmtry potential, we kept the isoscalar potential fixed 
while changing the symmetry potential from the momen- 
tum dependent symmetry potential to the momentum 
independent symmetry potential and keep the symmetry 
energy fixed |22j. Fig. [1] shows the momentum dependent 
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dpi/dt = -V r U(n) + qiE, 

dn/dt = p. l /^m 2 +p 2 + \7 p U(r l , Pt ), (2) 

where qi is the charge of particle, E is the Coulomb field 
of particle felt. The incomprcssibility of symmetric nu- 
clear matter at normal density is set to be 211 MeV. Ac- 
cording to essentially all microscopic model calculations, 
the EOS for isospin asymmetric nuclear matter can be 
expressed as 



E{p, 6) = E(p, 0) + E sym {p)6 2 + 0(S 4 ), 



(3) 



where E(p, 0) is the energy per nucleon of symmetric nu- 
clear matter, and E sym (p) is the nuclear symmetry en- 
ergy. With the single particle potential U(p,6,p,r), for 
a given value x, one can readily calculate the symmetry 
energy E sym (p) as a function of density. Because the pur- 
pose of present studies is just to see how large the effect 



FIG. 1: (Color online) The momentum dependent symmetry 
potential (MD) as a function of momentum and the momen- 
tum independent symmetry potential (MID) as a function of 
density (upper x-coordinate) from the Gogny interaction [22] . 



symmetry potential (MD) and the momentum indepen- 
dent symmetry potential (MID) from the Gogny inter- 
action. We can see that at higher densities strength of 
the momentum independent symmetry potential (MID) 
is always larger than that of the momentum dependent 
symmetry potential. In the following, we give our results 
of the momentum dependence of nuclear symmetry po- 
tential on nuclear transverse and elliptic flows at higher 
densities (The density reached in the 132 Sn+ 124 Sn at a 
beam energy of 400 MeV/nucleon is about 2 times satu- 
ration density |23j). 
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FIG. 2: (Color online) Rapidity distribution of free neutron 
to proton ratio n/p in the reaction 132 Sn+ 124 Sn at a beam 
energy of 400 MeV/nucleon and an impact parameter of 5 fm 
with and without momentum dependence of nuclear symme- 
try potential, signed with MD and MID, respectively. 



FIG. 3: (Color online) Neutron and proton transverse flows 
analysis in the reaction Sn+ Sn at a beam energy of 400 
MeV/nucleon and an impact parameter of 5 fm with and with- 
out momentum dependence of nuclear symmetry potential, 
signed with MD and MID, respectively. 



III. RESULTS AND DISCUSSIONS 

We first study the neutron to proton ratio n/p of free 
nucleons as a function of rapidity as shown in Fig. [2J It 
is seen the case with momentum dependence of nuclear 
symmetry potential causes lower neutron to proton ratio, 
whereas the case without momentum dependence of nu- 
clear symmetry potential causes higher neutron to proton 
ratio, especially for nucleons at large rapidities. Because 
the momentum dependence of nuclear symmetry poten- 
tial decreases the strength of nuclear symmetry potential 
at high densities or high nuclconic momenta as shown in 
Fig. [TJ The small symm etry potential decreases the free 
neutron to proton ratio [22|. From this plot, we can see 
that while using the neutron to proton ratio n/p of free 
nucleons to probe the symmetry energy, one should keep 
in mind that this observable is sensitive to the the mo- 
mentum dependence of nuclear symmetry potential used, 
and thus may cause uncertainties while compared with 
the experimental data. 

Fig. [3] shows neutron and proton transverse flows anal- 
ysis in the reaction 132 Sn+ 124 Sn at a beam energy of 
400 MeV/nucleon and an impact parameter of 5 fm with 
and without momentum dependence of nuclear symme- 
try potential. We can see that without (with) momentum 
dependence of nuclear symmetry potential, the strength 
of the nucleoli flow (especially neutron flow) increases 
(decreases). This is understandable since the momen- 
tum independence of nuclear symmetry potential overall 
increases the strength of nuclear symmetry potential as 
shown in Fig. [I] thus neutrons are repelled more strongly 
than the case with the momentum dependent symmetry 



potential. From Fig. [3] we can also see that effect of the 
momentum dependence of nuclear symmetry potential on 
the proton flow is less evident owing to the Coulomb po- 
tential added on protons. The Coulomb potential (always 
repulsive for protons) decreases the strength of symmetry 
potential (attractive for protons at high densities) added 
on protons. 
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FIG. 4: (Color online) Difference of proton and neutron trans- 
verse flows analysis in the reaction 132 Sn+ 124 Sn at a beam 
energy of 400 MeV/nucleon and an impact parameter of 5 fm 
with and without momentum dependence of nuclear symme- 
try potential, signed with MD and MID, respectively. 



The difference of proton and neutron transverse flows 
has been studied previously and was shown to be sensitive 
to the symmetry energy |4[ • In order to show if the differ- 
ence of proton and neutron transverse flows is sensitive to 
the momentum dependence of nuclear symmetry poten- 
tial, we plot Fig. |4j the difference of proton and neutron 
transverse flows analysis in the reaction 132 Sn+ 124 Sn at a 
beam energy of 400 McV/nuclcon and an impact param- 
eter of 5 fm with and without momentum dependence of 
nuclear symmetry potential. From this plot, we see that 
the difference of proton and neutron transverse flows is 
insensitive to the momentum dependence of nuclear sym- 
metry potential. It is thus better to be used to probe the 
symmetry energy since the the momentum dependence 
of nuclear symmetry potential is still an open question 
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FIG. 5: (Color online) Neutron proton differential transverse 
flow analysis in the reaction 132 Sn+ 124 Sn at a beam energy of 
400 MeV/nucleon and an impact parameter of 5 fm with and 
without momentum dependence of nuclear symmetry poten- 
tial, signed with MD and MID, respectively. 



The neutron proton differential transverse flow also has 
been studied extensively and was shown to be sensitive to 
the symmetry energy [3|, l22j . To see if neutron proton dif- 
ferential transverse flow is sensitive to the momentum de- 
pendence of nuclear symmetry potential, we give Fig. [5J 
the neutron proton differential transverse flow analysis 
in the reaction 132 Sn+ 124 Sn at a beam energy of 400 
MeV/nucleon and an impact parameter of 5 fm with and 
without momentum dependence of nuclear symmetry po- 
tential. From this plot, we see that the neutron proton 
differential transverse flow is sensitive to the momentum 
dependence of nuclear symmetry potential as shown in 

The neutron-proton differential transverse flow is de- 
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where N(y), N n (y) and N p (y) are the number of free 
nucleons, neutrons and protons, respectively, at rapidity 
y; pf(y) is the transverse momentum of the free nucleon 
at rapidity y; Wi = 1 (—1) for neutrons (protons); and 
(Pn(y)) an d (Pp{y)) are respectively the average trans- 
verse momenta of neutrons and protons at rapidity y. 
One can see from Eq. (TJJ that the constructed neutron- 
proton differential transverse flow depends not only on 
the proton flow and neutron flow but also on their rela- 
tive multiplicities. Therefore the neutron-proton differ- 
ential transverse flow is not simply the difference of the 
neutron and proton transverse flows, it in fact depends 
also on the isospin fractionation at the rapidity y. If 
neutrons and protons have the same average transverse 
momentum in the reaction plane but different multiplici- 
ties in each rapidity bin, i.e., (p£(j/)) = {Pp{y)) = {p x {y)), 
and N n (y) ^ N p (y), then Eq. @ is reduced to 



M 



N n (y)-N p (y) 
N(y) 



{p x (y)) = S(y)-(p x (y)), (5) 



reflecting effects of the isospin fractionation. On the 
other hand, if neutrons and protons have the same mul- 
tiplicity but different average transverse momenta, i.e., 
N n (y) = N p (y) but (j%(y)) ± (p x p (y)), then Eq. © is 
reduced to 



1 



*£-p(v) = 5(<p£(v)> -<$(!/)»■ 



(6) 



In this case it reflects directly the difference of the neu- 
tron and proton transverse flows. Because the effect of 
the momentum dependence of nuclear symmetry poten- 
tial for neutron to proton ratio is very large (shown in 
Fig. [2]) and the difference of the proton and neutron 
transverse flows has no such combination, we see larger 
effect of the momentum dependence of nuclear symme- 
try potential on the neutron proton differential transverse 
flow than the difference of proton and neutron transverse 
flows. It is thus better to probe the symmetry energy by 
using the difference of neutron and proton flows since the 
momentum dependence of nuclear symmetry potential is 
still an open question [12| ■ 

The elliptic flow V2(y,pt), which is derived as the sec- 
ond coefficient from a Fourier expansion of the azimuthal 
distribution N(<p,y,p t ) = v (l + v\cos{4>) 4- 2u 2 cos(20)), 
can be expressed as 



v 2 =< 



Px - PI 



PI 



>, 



(7) 



where pt 
I2II2F 



pI 



p 2 is the transverse momentum 
The difference of proton and neutron elliptic 



IV. CONCLUSIONS 
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FIG. 6: (Color online) Difference of proton elliptic flow and 
neutron elliptic flow analysis in the reaction 132 Sn+ 124 Sn at 
a beam energy of 400 MeV/nucleon and an impact parameter 
of 5 fm with and without momentum dependence of nuclear 
symmetry potential, signed with MD and MID, respectively. 



Based on the IBUU04 transport model, effect of the 
momentum dependence of nuclear symmetry potential 
on nuclear transverse and elliptic flows in the neutron- 
rich reaction 132 Sn+ 124 Sn at a beam energy of 400 
MeV/nucleon is studied. It is found that the momen- 
tum dependence of nuclear symmetry potential affects 
the rapidity distribution of the free neutron to proton 
ratio, neutron flow and proton flow as a function of rapid- 
ity. The momentum dependence of nuclear symmetry po- 
tential affects neutron-proton differential transverse flow 
more evidently than the difference of neutron and pro- 
ton transverse flows as well as the difference of proton 
elliptic flow and neutron elliptic flow. Therefore it is 
better to probe the symmetry energy by using the differ- 
ence of neutron flow and proton flow since the momen- 
tum dependence of nuclear symmetry potential is still an 
open question. And it is better to probe the momen- 
tum dependence of nuclear symmetry potential by using 
the neutron-proton differential transverse flow and the 
rapidity distribution of the free neutron to proton ratio. 



flows is also shown to be sensitive to the symmetry energy 
[251 |26| . It is thus also interesting to see if the difference 
of proton elliptic flow and neutron elliptic flow is sensi- 
tive to the momentum dependence of nuclear symmetry 
potential. Fig. [6] shows the difference of proton ellip- 
tic flow and neutron elliptic flow analysis in the reaction 
132 Sn+ 124 Sn at a beam energy of 400 MeV/nucleon and 
an impact parameter of 5 fm with and without momen- 
tum dependence of nuclear symmetry potential. Again, 
we see that the difference of proton elliptic flow and neu- 
tron elliptic flow is not sensitive to the symmetry energy. 
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